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The rheological response under simple shear of an active suspension of Escherichia coli is de¬ 
termined in a large range of shear rates and concentrations. The effective viscosity and the time 
scales characterizing the bacterial organization under shear are obtained. In the dilute regime, we 
bring evidences for a low shear Newtonian plateau characterized by a shear viscosity decreasing 
with concentration. In the semi-dilute regime, for particularly active bacteria, the suspension dis¬ 
play a ’’super-fluid” like transition where the viscous resistance to shear vanishes, thus showing that 
macroscopically, the activity of pusher swimmers organized by shear, is able to fully overcome the 
dissipative effects due to viscous loss. 


Owing to its relevance in medicine, ecology and its im¬ 
portance for technological applications, the hydrodynam¬ 
ics of active suspensions is at the centre of many recent 
fundamental studies [1, 2]. In nature, wide classes of liv¬ 
ing micro-organisms move autonomously in fluids at very 
low Reynold numbers [3]. Their motility stems from a va¬ 
riety of propulsive flagellar systems powered by nanomo¬ 
tors. For bacteria such a bacillus subtilis or E.coli, the 
propulsion comes from the rotation of helix shaped flag¬ 
ella creating a propulsive force at the rear of the cell body 

[4] . Consequently, many original fluid properties stem 
from the swimming activity [5-11]. Due to hydrodynamic 
interactions bacteria may produce mesoscopic patterns of 
collective motion sometimes called “bio-turbulence” [IS¬ 
IS]. In a flow, these bacteria may organize spatially [12] 
and under shear, for pusher-swimmers, the swimming ac¬ 
tivity yields the possibility to decrease the macroscopic 
viscosity to values below the suspending fluid viscosity 

[5] . In the dilute regime, kinetic theories via a sim¬ 
ple account of the dominant long range hydrodynamic 
held [19-21] provide closed-forms for shear viscosity as 
a function of shear rate. Remarkably, at low shear rate, 
these theories predict a Newtonian plateau with a vis¬ 
cosity decreasing linearly with concentration [19-21]. On 
the other hand, phenomenological theories were also pro¬ 
posed to describe macroscopically active suspensions via 
a coupling of hydrodynamic equations with polar and/or 
nematic order parameters ]2, 5, 6, 22-24]. A striking 
outcome of these theories is that for a set of coupling pa¬ 
rameters rendering essentially a high swimming activity, 
a self-organized motive macroscopic flow may show up 
in response to shear ]22-24]. This onset of a dissipation¬ 
less current is described in analogy with the super-fluidity 
transition [22, 23] of liquids. Experimental evidences for 
viscosity reduction to values below the suspending fluid 
viscosity were brought for Bacillus subtilis [8] and E. coli 
[25] suspensions. However, no full rheological characteri¬ 
zation (i.e. viscosity versus shear rate) under steady and 


uniform shear exists. Moreover, these pioneering experi¬ 
ments did not provide evidence for the low shear viscous 
plateau which is at the core of all theoretical predictions 
in the dilute regime. Finally, the phenomenological pre¬ 
dictions for the non-linear regime have remained so far 
unobserved. Noticeably, for unicellular algea, viewed as 
’’puller” swimmers, the predicted low shear rate increase 
of viscosity was measured experimentally [9]. In this let¬ 
ter, in addition to a full rheological characterization of 
an E. coli suspension, we provide, in the dilute regime 
and at low shear rate, experimental evidences for a linear 
decrease of the apparent viscosity with bacteria concen¬ 
tration. We also explore regimes of higher concentration 
and describe the conditions where we observe a transi¬ 
tion to a dissipation-free macroscopic flow. 

The active fluids considered here, are prepared out of 
two strains of wild type E. coli (ATCC9637 and RP437) 
suspended into a minimal medium where the bacteria 
are still motile but do not divide. ATCC9637 is cultured 
overnight at 25 °C in LB medium shaken at 240 rpm. 
RP437 is cultured overnight at 30 °C and shaken at 240 
rpm in M9 minimal medium supplemented with 1 mg/ml 
casamino acids and 4 mg/ml glucose. Next, the culture 
is washed twice by centrifugation (2300 g for 10 min) 
and the cells are re-suspended into a motility medium 
containing 10 mM potassium phosphate pH 7.0, 0.1 mM 
K-EDTA, 34 mM K-acetate, 20 mM sodium-lactate and 
0.005 % polyvinylpyrrolidone (PVP-40). To avoid bacte¬ 
rial sedimentation, the suspension is mixed with Percoll 
(1 vol/1 vol). The bacteria concentration n is represented 
by its volume fraction (j) = n/Vb, where 14 is the bacteria 
body volume chosen as the classical value 14 = Ipm/. 
Shear stress is measured in a low-shear couette rheome¬ 
ter (Contraves 30) designed especially for probing low- 
viscosity fluids. The inner bob (radius Ri = 5.5 mm, 
length 8 mm and underside cone angle 20° ) is suspended 
by a torsion wire into a cup (inner radius Ri = 6 mm). 
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FIG. 1. Shear stress response for an E.coli suspension 
(ATCC9637 strain, T=25°C'). (A) Shear stress E rescaled 
by the applied shear rate 7 during the rotation to display an 
effective viscosity E/7 in the sheared regime. Grey and black 
lines: fluid without bacteria (</ = 0). Colored solid lines: 
fluid with bacteria (</ = 0.67%). Various 7 are applied rang¬ 
ing from 64 s“'^ (dark blue line) down to 0.022 s“^ (dark red 
line). (B) Relative viscosity rj/rjo averaged over three realisa¬ 
tions as a function of 7 (□: (f> = 0.11%, O- 4> = 0.21%, A: 
cp = 0.44% and (): (j> = 0.67%). The solid line is an adjustment 
by the Carreau law: rf/rio = 1.08 — 0.795/(1 4- (7/0.6)^). The 
vertical dashed line shows 7^, below this shear rate rj{-y)/rio 
is less than 1. (C) Values of the plateau viscosity rjp/rio as 
function of the bacteria volume fraction (p for very low shear 
rates (Q: 7 = 0.022 s"\ □: 7 = 0.04 s"\ A: 7 = 0.075 s"^). 


The cup rotates at an angular rate controlled by a com¬ 
puter. The corresponding shear rate is 7 = ■ The 

central feature, making this instrument very precise for 
low stress measurements, is that the central bob is kept 
fixed by a feedback counter rotation of the suspending 
wire. The instrument measures the compensating torque 
required to keep the torsion wire at its null position. The 
torque is then converted into shear stress every 0.7 s. 
Importantly, due to the small surface area between the 
fluid and the air, the flux of O 2 is insufficient to compen¬ 
sate the amount of O 2 consumed by the bacterial activity. 
To avoid bacteria suffocation and consequently, a severe 
drop of activity, we supplement the suspension with L- 
Serine, an amino-acid allowing the bacteria to keep a sig¬ 
nificant swimming activity in absence of oxygen [27, 28]. 
Therefore, in the early instants of the measurements, the 
bacteria are still in oxygenated conditions; but since they 
consume the oxygen, their mean velocity and diffusion 
coefficient decrease and stabilize within about 10 min. 
Consequently, we observe a continuous increase in the 
suspension viscosity until a constant value is reached. 
Then, by metabolizing L-serine bacteria sustain a con¬ 
stant activity lasting for few hours (see Supporting Ma¬ 
terial (SM)). To obtain a full rheogram as displayed in 
Fig.IB, the following protocol is used. A volume (1.25 
ml) of the suspension is poured into the rheometer’s cup 
and then the bob is set into place. After 30 s of rest, 
the cup is rotated for 30 s at a steady state shear rate. 
The rotation is then stopped for 30 s. These steps are re¬ 
peated with increasing shear rate values. Once the high¬ 
est 7 is reached, the procedure is repeated in decreasing 
order to verify the reversibility of the viscous response. 

Rheology measurements were first performed to ob¬ 
tain viscosity at volume fractions between (p = 0 . 1 % and 
0.67%. In Fig. lA, we display the stress responses ob¬ 
tained for a suspension at a given (p, for various shear 
rates. For the suspending fluid alone (a Newtonian fluid, 
r]Q = l.dmPa.s), stress-time responses, at the start or at 
the stop of the applied shear, are fast and correspond 
to the device compliance (grey and black lines). A sim¬ 
ilar behaviour is observed for the suspensions probed at 
high shear rates which moreover display a viscosity higher 
than r]o as observed classically for suspensions of passive 
particles. However, at low shear rates, a strikingly differ¬ 
ent behaviour is observed. When shear starts, the stress 
jumps to the value measured in the absence of bacteria 
and then, after an exponential decrease, lasting for a few 
seconds, a steady effective viscosity - 77 - is reached. When 
shear stops, the stress decreases abruptly and eventually 
changes sign. Finally, the stress relaxes exponentially to 
0 with a characteristic time t~ not very different from 
the time scale, t+, needed to reach a steady viscous re¬ 
sponse under shear (data on Fig.2C). In this last stage, 
the bacterial motion induces a motive stress on the inner 
bob. 

Figure IB shows the suspension viscosity 77 as a func- 
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tion of 7 for different volume fractions ranging from 
(j) — 0.11 % (1.1 X 10® bact/mL) up to ^ = 0.67 % 
(6.7 X 10® bact/mL). We observe the three regimes pre¬ 
dicted by the theories [19, 20]; (i) at high shear rates 
(7 > ls“^), the active contribution to viscosity is negligi¬ 
ble and a Newtonian plateau appears akin to suspensions 
of passive particles of the same shape; (ii) below a criti¬ 
cal shear rate value jc < 1-5 s“^, the suspension viscosity 
is lower than the suspending fluid viscosity; (iii) at low 
shear rates (7 < 0.1 s“^), an ’’active viscous plateau” 
rip{(f>) appears. Futhermore, the theories also predict a 
linear dependance of rjp with (j) given by: 

^ = l + iF(f)</ ( 1 ) 

Vo tc 

where K cx (A—B^) ; tc is the time taken by a bacterium 
to drag the fluid over its size and t characterizes the di¬ 
rectional persistence of a swimming trajectory [26]; A 
and B depend solely on the bacterium shape [19-21]. As 
shown in Fig.1C, rjpi'P) decreases linearly with the con¬ 
centration as ^ = \+K(j) with K ~ —120±10, as long as 
the shear rate is sufficiently low (i.e. in the range 0.022 
- 0.075 s“^). The experimental results are very consis¬ 
tent with active suspension theories for pusher swimmers 
in the dilute regime. Noticeably, within the framework 
of closed form theories established in the dilute regime, 
such a complete determination of the viscous response 
can be used to give an estimation of the microscopic bac¬ 
terial activity (See SM for an explicit derivation of the 
dipolar strength and other microcopic parameters via the 
Saintillan [20] kinetic model). 

We next compare the viscous response of the bacteria 
for two different activity levels. When both O 2 and 
L-serine are present in the suspension, one obtains 
an ’’hyper-activated” regime characterized by high 
values of bacteria diffusivity and swimming velocity 
(D = 7 10-iimVs and Vb = 2Sv.m/s for ATCC9637). 
After 10 min, when O 2 is consumed, the motility is 
maintained by the metabolizaton of L-serine. In this 
case one has a lower activity with D = 1.2 10“^^m®/s 
and Vb = 20fim/s (See SM Fig.2). The first regime, 
only lasts few minutes and the measurement were thus 
restricted to a single 7 sufficiently low (here 0.04s“^) 
to estimate ryp, the active plateau viscosity. In the 
’’hyper-active” domain, r]p is again found for (j) < 0 . 6 % 
to decrease linearly with but with a larger slope 
K' ~ —200 ± 3 (see Fig. 2B). This results demonstrates 
that the slope K is strongly related to the bacterial 
activity. Next, these experiments were repeated with 
the second E-coli strain (RP437). A linear rjp vs (j) 
relation is also found but with a larger negative slope 
K” ^ —259± 13 (red squares in Fig. 2B) demonstrating 
a correlation between the bacterial characteristics and 
the rheological response at low shear rate. 
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FIG. 2. (A) Shear stress response E rescaled by the shear rate 
7 for the ATCC9637 strain (blue lines) and RP437 strain (red 
lines). All experiments are performed with 7=0.04 s“^. In 
some cases, the stress reponse can reach negative values. (B) 
Variation of the viscosity — as function of the volume frac- 
tion of bacteria (f) in oxygenated conditions (filled symbols) 
and deoxygenated conditions (empty symbols). Dashed lines 
are meant as guides only. (C) Relaxation time Tr obtained by 
adjusting exponentially the stress relaxation at the start of 
the shear (rotation) and at the end of the shear (relaxation) 
as function of the mean distance (^ = s) between bacte¬ 
ria (empty symbols: deoxygenated conditions, filled symbols: 
oxygenated conditions, blue symbols: ATCC9637 strain, red 
symbols: RP437 strain). 


We finally increase the number of bacteria in the solu- 
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tion to values above 0.6%. For the ATCC9637 strain in a 
medium that does not contain oxygen (empty symbols in 
Fig. 2B), we observe that the viscosity becomes constant 
and independant of (p for p > 0.7% with rjp/riQ ^ 0.2. In 
oxygenated conditions and for highly motile batches, we 
measure that the viscosity ijp/rjo also reaches a constant 
value beyond </> ~ 0.5% and p ~ 0.4% for ATCC9637 
and RP437 respectively. For the three cases, the ac¬ 
tive plateau viscosity becomes thus independent of the 
concentration over a significant domain (from tp = 0.6% 
up to 2.4%). In addition under ’’hyper-activated” condi¬ 
tions, one observes a viscous response reaching zero (see 
Fig. 2A), meaning that the local viscous dissipation is 
macroscopically entirely compensated by the swimming 
activity. Moreover, with the very active RP strain, nega¬ 
tive values for rjp (see Fig. 2A and 2B) could be obtained 
at the edge of the transition. These results point towards 
the idea of an organization process triggered by the shear 
flow which may last even when the shear ceases. In Fig. 
2C, we plot the relaxation times and t~ as function 
of the mean distance ^(= (p~^^^) between bacteria. The 
transition between the dilute and the semi dilute regimes 
occurs here for a distance ^ ~ 6/im. In fact, below this 
length, the relaxation time is found to decrease linearly 
with ^ (while the viscosity rjp remains constant). The 
result seems to be consistent with the phenomenological 
picture proposed for active polar swimmer and with the 
clustering interaction recently observed in absence of flow 
[17], and suggests that collective effect becomes impor¬ 
tant for distances shorter than 6/im. 

In conclusion our experiments show that the bacterial 
activity has a measurable influence on shear viscosity. 
We brought direct experimental evidence, at low shear 
rates, for an active viscous plateau which value decreases 
linearly with concentration (for (p < 0.3%). This con¬ 
firms a central prediction for active pusher suspensions 
in the linear kinetic regime. The most striking feature 
of the rheological response is indeed the emergence of a 
viscous-less ’’super-fluidity” regime (?/ ~ 0 or even lower). 
Presently, there is no ah initio micro-hydrodynamic cal¬ 
culation describing the impact of bacteria interactions 
and possibly, the influence of collective organization on 
the macroscopic rheology. Therefore, one has to rely on 
phenomenological arguments to identify the possibly rel¬ 
evant macroscopic hydrodynamic contributions [2, 22]. 
In this framework, several authors [22-24] investigated 
the responses to shear of active polar particles and re¬ 
markably, they predicted the possibility of a transition to 
a ” zero-viscosity” regime when the activity is increased. 
This result can also be cast in the framework of recent 
experimental works pointing out the possible use of bac¬ 
terial motion to drive mechanical devices [33-35]. Actu¬ 
ally this report goes in this direction as we show, at least 
in principle, that rotational macroscopic power could be 
extracted from the swimming activity as for a rotatory 
motor [24]. Finally, such a strong viscosity reduction 


may be a crucial element when considering macroscopic 
transport and particle dispersion in porous systems or in 
capillary networks, a central question to many applica¬ 
tions involving bacterial fluids. 
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SUPPLEMENTAL MATERIAL OF THE MANUSCIPT : 
"TURNING BACTERIA SUSPENSIONS INTO A 
"SUPERFLUID "" 

Variation of the viscosity rj from oxygenated to deoxygenated 
conditions 



Figure 1. Variation of the viscosity rj from oxygenated to deoxy¬ 
genated conditions at low shear 7 = 0.04 The ratio of visco¬ 
sities ^ of a suspension containing a volume fraction 0 = 0.21 % 
of bacteria (ATCC9637) is initially of the order of ~ 0.6 and reaches 
a constant value of ~ 0.76 after 20 min. This increase is due to the 
consumption by the bacteria of the O 2 initially dissolved in the solu¬ 
tion. After 20 min, the remaining activity is due to the metabolization 
of L-serine [1]. Inset : To demonstrate the relation between visco¬ 
sity variation and O 2 consumption, the same measurement is perfor¬ 
med with a solution containing twice the concentration of bacteria 
at (f> = 0.44 %. The data collapse observed when the time t is mul¬ 
tiplied by the concentration 0 demonstrates that it takes half of the 
time to reach the stabilized viscosity plateau suggesting a constant 
O 2 consumption rate. Thus, our data allow an estimation for the oxy¬ 
gen consumption rate by the bacteria. In this buffer we find a rate of 
« 2.5 X 10"* dioxygen molecules/s/bacterium which is close to the 
value reported by Douarche and coworkers [1]. 


Evolution of the microscopic characteristics of bacteria motility 



Figure 2. Evolution of the microscopic characteristics of bacte¬ 
ria motility. To confirm that the viscosity increase is due to a mo¬ 
dification of the bacteria swimming characteristic, a suspension of 
bacteria ((p = 0.3 %, ATCC9637) is either placed in a sealed sample 
(i.e. between two glass plates separated by a distance of 270/rm and 
sealed by a GeneFrame gaslight sealing system) or in a micro-fluidic 
channel made with an air permeable material. The bacteria are tra¬ 
cked at 125 fps and the mean square displacement is next com¬ 
puted as a function of time. Hence we extract the mean velocity 
Vo{t) and the diffusion coefficient D{t). In the micro-fluidic cell, 
where no shortage of O 2 occurs D and Vb are constant and close 
to 140 and 28 /im.s“^ respectively. In the sealed sample, 

the bacteria consume the O 2 until exhaustion and turn deoxygented 
in presence of L-serine. This corresponds to a constant drop of both 
D and Vo which finally plateau at values : D = 13 /rm^.s“^ and 
Vo = 20 The time characterizing the directional persistence 



Estimation of the dipolar strength \ao\, the thrust force /, and 
the effective size of a single bacterium I, via a slender body 
model. 


Here we describe a method to obtain from the rheological 
responses - the rheogram 77 ( 7 ) and the stress relaxation time 
(Tr) - the dipolar strength |(To| of the bacteria. The latter can 
be seen as the work produced by one bacterium to move the 
fluid over its own size at a swimming velocity Vq. To do so, 
we use the theoretical and numerical results of Saintillan [2] 
and we futhermore make the assumption that the relaxation 
time of the signal gives an estimation of the t the characte¬ 
ristic time of directional persistence of the trajectories. In this 
article, Saintillan computes the contribution of active slender 
rod-like particles on the total viscosity rj for a suspending 
fluid of viscosity 779 . First, we use jc the shear rate for 
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which starts to decrease below 1 ) to estimate tc, the 
typical time for a bacteria to drag the fluid along its length 1. 
Fig. 4 in Saintillan’s article [2] displays for different values 
of ^ the corresponding values of 7 cfc- In the actual range 

*c 

of operation of our experiment, we realize that a linear fit 
of the numerical data with ^ ~ 9 jctc allows to obtain a 
useful analytical relation between these variables. From the 
experimental measurements ~ 2.4 ± 0.3 s (See Fig. 2c) 
and 7c ~ 1-5 ± 0.25 s"! (See Fig. lb), we obtain a value 
fc = 0.42 ±0.05 s. 

Finally, using the Eq.(18) of the article of Saintillan [2] and 
considering the bacteria as slender rod-like particles (taking 
the Bretherton constant /3 = 1), we obtain the following ana¬ 
lytical expression for rjp : 


slopes K and K' corresponding to two different swimming 
velocities Vb and when bacteria swim with and without 
oxygen respectively. Since K = — X), where 

X = ^ and tc oc XjV then for the two bacteria velocities 

Vo and Vq, we have ^ . For the ATCC9637 strain, 

under oxygenated conditions we have Vq = 20 and 

X = ^ = 5.6 ± 0.6. When the cell activity is maintained 
by the metabolization of both O 2 and L-Serine (oxygenated 
condition), we measure a velocity Vq = 28 /im.s“^ and no 
significant change of is observed (See Fig. 2c). The ratio 
of the slopes is thus : ^ ~ 0.41 ± 0.08. This results is in 
agreement with our experimental measurement of the ratio ; 

=0.56 ±0.05. 


m 15 ? 7 o ^ tc 


where n = — is the number of bacteria per unit of volume 
and Vb is the volume of a single bacteria. Since the empirical 
parameter K is defined as ^ = 1 ± the dipolar strength 
(Jo can be obtained from the microscopic parameters r^, tc and 
K : 


kol 


I'oKvb'qo 
Ate - Tr 


( 2 ) 


Using the values K ~ —120 ± 10, Tr ~ 2.4 ± 0.3s and 
tc ~ 0.42 ± 0.05s measured for the ATCC9637 strain when 
supplemented by L-Serine and Vb = A y^m^, we obtain : 
|(To| = 3.8 ± 1.0 10“^® Pa.m® By using the expression 

.3 

given by Saintillan [2] : tc = in{Si) where I and d 

are the dipolar length and the diameter of the bacterium 
respectively, we find that the corresponding dipolar length 
I ~ 20.7 ± 3 /im corresponding to the flagella bundle length 
plus the E. coli body size [3, 4]. The dipole force / = \<Jo\/l 
can also be estimated. For our measurements, we find a value 
of / ~ 0.23 ± 0.04 pN. This agrees with recent experiments 
for swimming E. coli reporting flagella thrust forces / ^ 0.57 
pN [5], / ~ 0.41 pN [6] and / ~ 0.42 pN using PIV 
measurements [7]. 


From the previous slender-body model, one can also 
predict (everything else being constant) the ratio of two 


Other models can be considered to estimate the microsco¬ 
pic property of the bacteria from rheological measurements. 
Haines and coworkers [ 8 ] for instance considered the cell as a 
rigid prolate spheroid and Ryan et al. [9] as a point-like dipole. 
They both predict a Carreau-like variation of the viscosity of 
the form : 


ryo ^ I (l±(i)2) 


(3) 


where I is a typical length associated to the bacteria. Our ex¬ 
perimental observations do agree with this type of variation 
(see solid line in Fig. lb). 
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